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Abstract 

We derive simple formulas connecting the generalized Wigner func- 
tions for s-ordering with the density matrix, and vice-versa. These for- 
mulas proved very useful for quantum mechanical applications, as, for 
example, for connecting master equations with Fokker-Planck equa- 
tions, or for evaluating the quantum state from Monte Carlo simula- 
tions of Fokker-Planck equations, and finally for studying positivity 
of the generalized Wigner functions in the complex plane. 

RIASSUNTO: In questo lavoro deriviamo semplici formule che con- 
nettono direttamente le funzioni generalizzate di Wigner con la rispet- 
tiva matrice densita. Queste formule sono molto utili in applicazioni 
quantomeccaniche, come, ad esempio, nel connettere master equations 
con equazioni di Fokker-Plank, o per determinare lo stato quantistico 
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da simulazioni di equazioni di Fokker-Plank, o, infine, per determinare 
la positivita delle funzioni di Wigner generalizzate sul piano complesso. 
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Since the Wigner's pioneering work Pj, generalized phase-space tech- 
niques have proved very useful in various branches of physics ||. As a 
method for expressing the density operator in terms of c-number functions, 
the Wigner functions often lead to considerable simplification of the quantum 
equations of motion, as for example, transforming operator master equations 
into more amenable Fokker-Planck differential equations (see, for example, 
Ref. ||). By the Wigner function one can express quantum- mechanical 
expectation values in form of averages over the complex plane (the classi- 
cal phase-space), the Wigner function playing the role of a c-number quasi- 
probability distribution, which generally can also have negative values. More 
precisely, the original Wigner function allows to easily evaluate expectations 
of symmetrically ordered products of the field operators, corresponding to the 
Weyl's quantization procedure ||. However, with a slight change of the orig- 
inal definition, one defines generalized s-ordered Wigner function W s (a,a), 
as follows 

W s (a,a)= [ ^e aI - SA+ i |A|2 Tr[D(A)£] (1) 

where the integration is performed on the complex plane with measure 
cP\ = dReXdlmX, D(a) = e aa] - aa denotes the displacement operator, and 
a and eft ([a, a)] = 1) are the annihilation and creation operators of the field 
mode of interest. Then, using the Wigner function in Eq. ([!]) one can evalu- 
ate s-ordered expectation values of the field operators through the following 
relation 

Tr[:{a^) n a m : s g] = J d 2 aW s {a,a)a n a m . (2) 

It is easy to show that the particular cases s = —1,0, 1 lead to antinormal, 
symmetrical, and normal ordering, respectively, in which cases the general- 
ized Wigner function W s (a,a) historically was denoted with the following 
symbols and names 

!^Q(a,a) for s=-l "Q-function" 
W(a,a) for s=0 (usual Wigner function) (3) 
P{a,a) for s=l "P-function" 

For the normal (s = 1) and antinormal (s = —1) orderings, the following two 
simple relations between the generalized Wigner function and the density 
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matrix are well known 

Q(a,a) = (a\g\a) , (4) 
g = J d 2 a P(a,a) \a) (a\ , (5) 

where \a) denotes the customary coherent state \a) = D(a)\0), |0) being the 
vacuum state of the field. Among the three particular representations (|3|), 
it is also well known that the Q-function is positively definite and infinitely 
differentiable (it actually represents the probability distribution for ideal joint 
measurements of position and momentum of the harmonic oscillator: see, 
for example, Ref. Qj). On the other hand, the P-function is known to be 
possibly highly singular, and the only pure states for which it is positive are 
the coherent states [0] . Finally, the usual Wigner function has the remarkable 
property of providing the probability distributions of the quadratures of the 
field in form of marginal distributions, namely 

J dlmaW(ae^,ae- i(p ) = (Rea|£|Rea)^ , (6) 

where \x),p stands for the eigenstates of the field quadrature = 
(a^e 1 ^ + h.c.)/2 (any couple of conjugated quadratures X^, X^ +7r /2, with 
[X$, X^ +7r / 2 ] = i/^i are equivalent to the position and momentum of a har- 
monic oscillator). Usually, negative values of the Wigner function are viewed 
as signature of a nonclassical state (one of the more eloquent examples is 
given by the Schrodinger-cat states |§ whose Wigner function is character- 
ized by rapid oscillations around the origin of the complex plane). From Eq. 
([!]) one can see that all s-ordered Wigner functions are related to each other 
through the convolution relation 

W s (a,a) = fd 2 pWAP,P) ,, 2 _ exp(-^-|a-/3| 2 ) (7) 

Eq. (0) shows the positiveness of the generalized Wigner function for s < — 1, 
as a consequence of the positiveness of the Q-function. From a qualitative 
point of view, the maximum value of s keeping the generalized Wigner func- 
tions as positive can be considered as an indication of the classical nature of 
the physical state. 
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In this paper we present three equivalent trace forms that connect s- 
ordered Wigner functions with the density matrix. They are the following: 



W„(a,a) 



7r(l — S 

2 



-e-i^H 2 Tr 



7T(1 -S) 



e^ |Q|2 Tr 



S + l\ -^aa~ -r^-aat 

- e 1 - 8 ^e 1 - 8 

s — 1 

s + 1 
s - 1 



e !+ s ge 1 + s 



^M 2 



tt(1 - s) 



Tr 



VT^J vi - s 



, (9) 
, (10) 

• (11) 



Eqs. (PHTTJD can be compared with the Cahill-Glauber formula M 



WJa,a) 



tt(1 - s) 



Tr < : exp 



1 -s 



(a — er)(a — a) 



(12) 



where the colons denote the usual normal ordering; Eq. (|TT| ) represents a 
generalization of the formula M 



2 

W^c^d 7 ) = — Tr \gD(2a) exp(z7r(A 



(13) 



Vice versa, the density matrix can be recovered from the generalized 
Wigner functions using the following expression 

t=£-J*<M.W***» ( T ^at)(l^i)"exp (*-„) . (14) 
The proof of our statements requires the following identity 



|0)(0|e 



a<9 Q 



e H 2 +afA-aA-i|A| 2 = 



(15) 



a=a=0 

which is proved in the Appendix. Then, through the following steps: 

d 2 X „t 



W s (a,a) 



7T- 



3 aA-aA+§|A|' 



Tr[D(\)g] 



d 2 \ „v 



oA— aA 



7T- 



Tr 



e at ^|0)(0|e a9 ^' 



1 2 +/3A-/3A+i(^l)|A| 2 



/3=/3=0 
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7T(1 - S) 

2 

7T(1 - s) 



Tr 



e a ^\0)(0\e^g 



/3=/3=0 



e !- s 



H 2 Tr 



2-nrat / l+_g 

1 - S 



one proves Eq. (0). Continuing from the last result we have 



WJa,a) 



Tr 



tt(I-s) 

1 _|_ s \ 5 atfl /l — s\ 5 ata 



1 - s 
2 



tt(i - s; 



1 + s 

e -T^H 2 Tr 



1 4- Q \ 

J—aa^ I ____ \ r^~- •<"' 

1 - S 



e 1 - s 



1 + s 



1 - s 



Q 



1 + S \l ata 2 



1-s 



1 + s\ 2 
1 - s 



I at. 



2 -i^H 2 Tr 



7T(1 - S) 

which proves Eq. (0). Eq. (|T0|) is derived using the following identities 

_2a_t/S+l\ ^ a f S + L\ _JSL. t ^ a 

e i- 3 a e 1 - 8 = e ei-* , 



s-1 



S-1 



s-1 



e i- s - a e i+s a =e i-^e !+ s 



s-1 



e !+ s 



As a check, from Eqs. (^)-|Tl|) one can easily recover the usual definition of 
the Wigner function ([[]) for s = 0, and Eq. (f|) for the Q-function (s = —1), 
namely 



W-i(a, a) 



-e-H 2 Tr 

7T 

-Q(a,a) 

7T 



-0 + ) ata e aa ge aai 



-e" |Q|2 Tr \\0)(0\e aa ge aat 

7T L 



The inversion formula (fH] ) is obtained using Eq. ([IT]) and the following 
formula H 



O = / — Tr[OD{a)]DUa) 

J 7T 



(16) 
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that holds true for any Hilbert-Schmidt operator O, and hence for a (trace- 
class) density matrix. One has 



s + 1\ 2 
1 - s 

d 2 a 



4 at, 



late. 



71 



-Tr 



1-s 5 



S + 1\2 

1 - s 



1 - s 



Q 



1 - s 



W s (a,a)—^— — -e 1 ^ 1 1 D 1 



71 



2ft 



(17) 



Hence, 
2 



: 



1 + s 
2 

1 + s 



l + s 



i„.f, 



-i \ —a 1 a 2a + 

d 2 aW,(a,a)e-^\ 2 { "^=" 

l + s 



2ft 



iV / T^s 

sat, 



1 — s\ ^ ata 



l + s 



m=a (\ - S\5 ata 



1 + S 



and then the result follows easily. In particular, for s = one has the inverse 
of the Glauber formula 



q = 2 J d 2 aW(a,a)D{2a)(-f a 



(18) 



whereas for s = 1 one recovers the relation ([|) that defines the P-function. 

The trace form of Eqs.(|9|- |T0| - |iT|) can be used for an analysis of positivity 
of the Wigner function, usually a quite difficult task, as confirmed in Ref. 
| |10| . In particular, from Eq. ([5]) one can immediately see that for s < 1 
(namely, with the only exception of the P-function) the s- Wigner function 

2 — 2 t 

can become negative, because the operator e 1 - saa pe 1 - a °' a is positive-definite, 

/ -4-1 \ a ^ a 

whereas the preceding factor [jzi ) is negative for s < 1, and positivity is 
guaranteed only for products of positive operators. On the other hand, from 
Eq. (|TTD one can easily see that there is always a state (the eigenstate of a) a 
with odd eigenvalue) that makes the s- Wigner function at a = negative 
for s < 0. 

The representations (3, |10|) for the generalized Wigner functions also pro- 
vide the easiest way to derive differential representations for boson operators 
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acting on a density matrix. By defining, analogously to Eq. ([!]), the gener- 
alized Wigner symbol for any operator O, 

W s (a,a\6) = J ^e aJ -^ + ^ 2 TT[D(X)d] , (19) 
from Eqs. (|9| JT0|) one immediately derives the relations 

W s (a,a\ag) = e -^H 2 Iz_f «^ e ^ H *W.(a, a) , 

a + ^-d^j W s {a,a) (20) 



^ s («,«|at^) = e — H (-^-9 Q )e- — H V S 



a, a) 



a-^d a )W s (a,a) , (21) 



and analogous relations for right multiplication by the boson operator. More 
generally, one can write a differential representation for any super- operator — 
i. e. right or left multiplication by an operator O — namely 

W s (a,a\Og) = F s [6-]W s (a,a) , W s (a,a\eO) = F s [-6]W s (a,a) , (22) 

where O- and -O denote left and right multiplication by the operator O, 
respectively, and F s are differential forms functions of a, a, d a and da with 
the following properties 

F s [d l d r \ = F s [6r\F s [dr\ , (23) 

F S [6 1 6 2 } = F S [-6 2 }F S {-6 1 ]. (24) 

[F s [-6 1 ],F 1 [6 2 -]] = 0, (25) 

F s [-6]=F s p.} . (26) 

The functional forms of the basic super-operators are summarized in Table |l|. 
The representations of -cr and -a can be easily obtained from those of Gl- 
and al- using identities (p6|) . Eq. (|25| ) is just the obvious statement that 
"left multiplication commutes with right multiplication" (for a and this 
corresponds to the identity [d a + K,a, da + na] = 0). Then, the differential 
representation of higher-order super-operators is easily obtained from the 
composition rules (53) and (|24]). Using the differential representation for 
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Super-operator 


F s 


a- 


a + 


at- 


a ■ 


- -^<9 

2 Q 


•a 


a ■ 


2 ^ 


•at 


a + ±2^(9 Q 


a • a) 




a 


2 + if* ( 1 + a<9 a + a^) + ( V ) * & 3 


a) ■ a 




a 


2 - i±2 ( 1 + a<9 a + aS^) + (i±2 ) 2 3 Q « 


a)a- 




a 


2 + i[(l - s)a<3a - (1 + - (1 + s) - i(l - s 2 )<9 aS ] 


■a)a 




a 


2 + - s)ad Q - (1 + s)acW- (1 + s) - |(1 - s 2 )<9 aS ] 



Table 1: Differential Wigner representation of some super-operators 

Bose super-operators, one can convert master equations into (possibly high 
order) Fokker-Planck equations. For example, the master equation of the 
damped harmonic oscillator (damping coefficient 7 and thermal photons n) 

dtQ = ~l^J l + + go) a — 2a^a^) — ^n(aa^ g + gaa) — 2a) ga) , (27) 

can be converted into the equivalent Fokker-Planck equation for the s-ordered 
Wigner function 

d t W s (a,a) = | [d Q a + c\a + (2n + l- s)d aS ] W s {a, a) . (28) 
For solving Fokker-Planck equations, one can use very efficient Monte-Carlo 



Green-function simulation methods (see, for example, Ref. |11|), choosing the 



parameter s such that both the Wigner function and the diffusion coefficient 



remain positive during the evolution. Then, from the inversion Eq. (14) 
one can recover the matrix elements (n\g\m) of the operator g in form of 
Monte-Carlo integrals of Laguerre polynomials. 

In conclusion, we have presented simple trace formulas that connect the 
generalized Wigner functions with the density matrix, and vice-versa, and 
we have shown how they can be practically used for: i) studying positivity 
of the generalized Wigner functions; ii) connecting master equations with 
Fokker-Planck equations; in) evaluating the quantum state in Monte Carlo 
simulations of Fokker-Planck equations. 
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Appendix 



Proof of identity fllS] 
From the relation 



gngm 
a a 



e H2 = 5 

c — u r , 



,ni 



(29) 



one has 



e at ^|0)(0|e a9 ' 

oo 

= E(« t r|0)(0|a"^ = EH^I = i - 



o=a=0 



M 1 



n=0 



e>- = E (« t ) n |0)(0| 

a=a=0 n,m=0 

oo 

n! 



On am 



nimi 



>l 5 



a=a=0 



n=0 



Hence, using the identities 



one obtains 



|0)(0|e 



e |a| 2 +aA-aA-i|A| 2 



a=a=0 

e aU e Q 



e -i|A| 2 e aA-aA e atA e at^| O ^ O | e a0 a 



e M 2 e -Aa 



o=o=0 



e -||A| 2 e atA e -aA = ^ 



(30) 



(31) 



(32) 
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